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o . 

f"^ We present a diagonalization method for generic matrix valued Hamiltonians based on a formal 

, expansion in power of ft. Considering ft as a running parameter, a differential equation connecting 

two diagonalization processes for two very close values of ft is derived. The integration of this 
differential equation allows the recursive determination of the series expansion in powers of ft for 
the diagonalized Hamiltonian. This approach results in effective Hamiltonians with Berry phase 
■ corrections of higher order in ft, and deepens previous works on the semiclassical diagonalization 

' of quantum Hamiltonians which led notably to the discovery of the intrinsic spin Hall effect. As 

physical applications we consider spinning massless particles in isotropic inhomogeneous media and 
show that both the energy and the velocity get quantum corrections of order ft 2 . We also derive 
formally to all order in ft the energy spectrum and the equations of motion of Bloch electrons in an 
external electric field. 
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I. INTRODUCTION 



Recently, a lot of works have shown the relevance of the Berry phase in semiclassical Physics. For instance, a new 
set of semiclassical equations of motion including a Berry phase correction was derived within a Lagrangian formalism 
to account for the adiabatic wave-packet evolution of an electron in an electromagnetic field (see also 0). A similar 
semiclassical Lagrangian approach for a light wave-packet in a inhomogeneous medium predicted an optical Hall effect 
||. However this Lagrangian formalism, although adapted to the evolution of a wave-packet, is only semi-classical by 
construction so that higher order quantum corrections are out of reach. Moreover, it was found that the subsequent 
derivation of the Hamiltonian presents some difficulties due to the presence of Berry- phase terms [H . 

The second difficulty has been solved by a diagonalization procedure with accuracy h for a generic matrix valued 
Hamiltonian 5] . This method results in an effective semiclassical diagonal Hamiltonian with Berry phase corrections 
■ as well as noncommutative covariant coordinates and momentum operators. The resulting generic equations of motion 
are also corrected by Berry phase terms. In the case of a Bloch electron in electromagnetic field 5] 6], this approach 
t-H ■ leads to the same equations of motion and the same magnetization as in [l[ . Therefore, the semiclassical diagonalization 
t-H ' of the quantum Hamiltonian not only shows that there is no trouble with the Hamiltonian in the presence of Berry 
phase terms, but also presents a practical short-cut to derive an effective Hamiltonian, noncommutative coordinates, 
momentum operators and equations of motion with Berry phase corrections. Moreover, this formalism, which in our 
opinion, better reflects the physical origin of the phenomena under consideration, is a general one that has been and 
could still be applied to several other systems in condensed matter or particle physics. 

For instance, for spinning particles, this approach led to the discovery of the intrinsic spin Hall effect or topological 
spin transport. Indeed it was already known that Berry phases, which are in fact spin-orbit couplings, modify 
semiclassical dynamics of spinning particles in electric Q and magnetic field [H as well as in semiconductor [§]. In 
this last case, the Berry phase in momentum space might be responsible for a transverse dissipationless spin-current 
in the presence of an electric field. In addition to that, spin-orbit contributions to the propagation of light in isotropic 
inhomogeneous media has been the focus of several other works Q [13] and has led to a generalization of geometric 
optics called geometric spinoptics flTlj j . 

The diagonalization method of Q allows to extend these studies to spinning particles interacting with static grav- 
itational fields and goes beyond previous approaches 12j HI- It led to the discovery of new couplings between the 



spin and magneto-torsion fields which could reveal an hypothetical torsion of space [141 ] [151 ] . 

Therefore, the diagonalization procedure of general quantum Hamiltonian has the advantage, compared to the 
semiclassical wave-packet Lagrangian approach, to gather several different systems into one general scheme. We find 
particularly interesting that our approach also gives to photons and electrons an equal footing. Indeed the spin Hall 
effect which offers promising applications both in spintronics and in spinoptics, has its origin, within our formalism, in 
the noncommutativity of the coordinates for both particles. From this point of view it is always legitimate to wonder 
whether an electronic phenomenon has its photonic counterpart. 

However, the diagonalization procedure mentioned above, despite its advantages, is limited to the semi classical 
level. It is thus legitimate to wonder about the possibility to go one step further, i.e. to find higher order quantum 
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corrections. The goal of the present article is to show how to go beyond the semiclassical approximation, a task 
beyond the reach of the Lagrangian formalism mentioned previously. The method we develop here is entirely new 
and different from the semiclassical diagonalization of [||, although based on it. It results in an effective diagonal 
Hamiltonian with Berry phase corrections as well as noncommutative coordinates and momentum operators (in the 
adiabatic approximation) all written as series expansions in h. 

More precisely, we propose a procedure to diagonalize recursively in power series of h a generic matrix valued 
Hamiltonian. The key element of this method relies on considering h as a running parameter. It allows to derive a 
differential equation connecting two diagonalization processes for two very close values of h.The recursive diagonal- 
ization can then be obtained by integrating this differential equation order by order in h. While this method does 
not allow an exact derivation of the diagonalized Hamiltonian, it has nevertheless the advantage to make formally 
possible the computation of quantum corrections of order higher than h (at least second order). To show this point 
and explain practically our technique, we study as a first example, a photon in an isotropic inhomogcncous medium 
at the second order in h, giving thus corrections to the semi classical equation we derived in [l4l ]. An important 
consequence of this computation is that the light velocity and the photon energy get a quantum correction of order 
h 2 . As a second different physical application, we also provide, formally to all order in h, the energy spectrum and the 
equations of motion for Bloch electrons in a constant external electric field. We show that the equations of motion 
are formally the same to all order in h, and that only the Berry connections get contributions in this expansion. 

The paper is organized as follows. In the next section we develop our formalism in the case of a generic matrix 
valued Hamiltonian. We derive the recursive differential equation in h for the diagonalized Hamiltonian. In section 3 
we discuss the semiclassical approximation and derive the generic equations of motion with Berry phases corrections. 
Section 4 is devoted to exemplify the method beyond the semiclassical level for physically relevant systems. Section 
5 is for the conclusion. 



II. RECURSIVE DIAGONALIZATION OF QUANTUM HAMILTONIAN 

In this section we consider a quantum mechanical system whose state space is a tensor product L 2 (7l 3 ) V with 
V some internal space. In other words, the Hamiltonian of this system can be written as a matrix Hq (P, R) of size 
dim V whose elements are operators depending on a couple of canonical variables P and R, the archetype example 
beeing the Dirac Hamiltonian with V = C 4 . Our goal is to describe a diagonalization process for this matrix valued 
quantum Hamiltonian Ho (P, R) recursively as a series expansion in powers of h. This expansion in powers of h gives 
the quantum corrections to the diagonalized Hamiltonian with respect to the classical situation h = 0. For example, 
the first order correction in h corresponds to the semiclassical approximation. Let us stress at this point that, by 
diagonalization, we always mean here a unitary transformation setting the Hamiltonian in a diagonal matrix form, 
the diagonal elements beeing operators depending on P and R. That is, we do not aim at finding the eigenvalues, 
but rather to derive the diagonal representation of Hamiltonians, that are usually relevant for the dynamics. 

We will derive the h expansion recursively in the following way. The Planck constant h is formally promoted to a 
dynamical parameter a in order to establish a differential equation connecting the two diagonalized Hamiltonians at 
h = a and h — a + da. The integration of this differential equation allows then the recursive determination of the 
different terms in the expansion of the diagonalized Hamiltonian in powers of a. 

A. Set up 

To start with, consider a quantum mechanical system with canonical variables P and R, where P is the generator 
of translations. It could be the usual momentum, or the magnetic translation operator of electrons in magnetic Bloch 
bands in solid state physics as discussed in the section devoted to the physical applications (see also @). 

Let consider the commutator between the canonical variables as a dynamical parameter a, that is 

[P,R] = -ia (I) 

Consider now the three following assumptions : 

1. The Hamiltonian of the system can be written as a matrix of a certain size Ho (P, R), that is a matrix whose 
coefficients depend on P and R. The typical example is the free Dirac Hamiltonian (a 4 x 4 matrix), but as shown 
in Q an electron in a periodic potential also fits in this set up. 

2. Assume moreover, that for each value of a, Hq (P, R) is exactly diagonalized through a matrix U a (P, R), i.e. 



U a (P, R) H Q (P, R) U+ (P, R) - e a (P, R) if [P, R] = 



—ta 



(2) 
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where e a (P,R) is a diagonal matrix. 

The index in L/ Q (P,R) underlines the fact that the diagonalization matrix depends a-priori on a since the diag- 
onalization process involves recombination of powers of P and R, and thus involves commutators of these variables 
(seed). 

3. The diagonalization is known when a — 0, that is when P and R commute. This assumption is crucial, since it 
allow to start a recursive diagonalization process in powers of a. This last assumption has also the advantage to be 
quite reasonable, since in many cases, the diagonalization is very easy to perform when [P, R] = 0. Actually, in this 
case the Hamiltonian can be seen as depending only on P for example, R being an external parameter. More about 
this point can be found in p|. 



B. The differential equation 

1. a as a parameter 

We will now consider a as a running parameter in order to find a relation between e a (P,R) and £ a +da (P>R)- 
But to facilitate the computation of e a (P, R) as explained in Q, we first need to write the Hamiltonian as well as all 
expressions involving the canonical operators R and P, in a symmetrized form. By this, we mean that whatever the 
prescription for the arrangement of the canonical variables (full symmetrization for instance or Weyl prescription) in 
the initial Hamiltonian, one will always rewrite it (and other expressions too) in a form where all the powers of P 
have been put half on the left and half of the right of the expression. Of course this symmetrization will introduce 
terms of order a due to the commutation relations. Note again that this is just a convention which facilitates the 
diagonalization. 

To start with, consider the diagonalization at the scale a 

U a (P, R) Hq (P, R) U+ (P, R) = e a (P, R) if [P, R] = -ia (3) 

and similarly for a + da. 

U a+da (P, R) H Q (P, R) U+ +dQ (P, R) = e a+da (P, R) if [P, R] =-i(a+ da) (4) 

Let us develop this last relation to the first order in da, 

(P, R) = U a H U+ + da (d a U a H U+ + U a H d a U+) 

where for clarity we have omitted the P and R dependence in the right hand side (note that here we assume the 
relation [P, R] = — i (a + da) and not [P, R] = — ia). This latter expression at the first order in da, reduces to 

e a+da (P,R) = U a H U+ + da(d a U a U+s a (P,R)+s a (P,R)U a d a U+) 

= U a H Q U+ + e a (P, R) (O a U a U+ + U a d a U+) da + [d a U a U+,e a (P, R)] da 

This differential equation is the main point of this paper. We now show how to compute explicitly each term of the 
right hand side, to make it tractable. 

2. Computation of d a U a (P, R) U+ (P, R) + U a (P, R) d a U+ (P, R) 

From the equality 

C/ Q+dQ (P,R)[/+ +dQ (P,R) =1 when [P,R] = -i (a + da) (5) 

one obtains 

da(a ct C/ ct (P,R)C/+(P,R) + C/ Q (P,R)9 Q C/+(P,R)) = 1 - C7„(P,R)C/+(P,R) (6) 
Note that here U a (P, R)C/+(P, R) j^l, since [P, R] = — i (a + da). As a consequence, we can ultimately rewrite : 

e a+da (P, R) = U a H U+ + e a (P, R) [l - U a U+] + [d a U a U+,e a (P, R)] da (7) 
we now need to rewrite U a HoU+ . 
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3. Computation of U a (P, R) Ho (P, R) (P, R) and final form for the differential equation. 

Here, the important point to keep in mind is that [P, R] = — i (a + da), so that 

U a (P,R)H (P,R)U+(P,R)^e a (P,R) (8) 

To compute U a H U+ we introduce two "fictitious" variables ,r, p, commuting with P and R and such that [p, r] = ida. 
We can consider r and p as an arbitrary couple of "small" canonical variables each of magnitude \J~da. As a 
consequence 

[P + p,R + r] = -ia (9) 

so that in virtu of Eq[2] we have the equality 

U a (P + p, R + r)H Q (P + p, R + r) C/+(P + p, R + r) = e a (P + p, R + r) (10) 

since P + p and R + r form a couple of canonical variables with commutator —ia. By the same trick, and for practical 
purpose, we can write : 

t/ Q (P + p, R + r)C/+ (P + p,R + r) =1 (11) 

As a consequence, using Eqs. 1101 and fTTI we can rewrite our differential equation Eq. [7] as : 

e a+da (P,R) = e Q (P + p,R + r) 

-U a (P + p, R + r)H (P + p, R + r) U+(P + p, R + r)-U a (P, R)H (P, R) U+(P, R) 
+e a (P, R) [U a (P + p, R + r)U+(P + p, R + r) - U a (P, R)U+(P, R)] 
+ [8 a U a (P, R) U+ (P, R) , e a (P, R)] da 

Now, we expand the R.H.S to the second order in powers of p and r. Since those two variables satisfy non trivial 
commutation relations, we choose to expand the R.H.S. on the basis Pi,Ti,PiPj,riVj ,PiTj for i ^ j, and P iri + 1 ~ i P i 
(Birkhoff-Witt factorization theorem). As a consequence, the term p^n or have to be rearranged in combination 

of ardpa ano j aazoa = y°. 

The expansion computation is similar to the one presented in Q (apart from some minor additional terms), and 
thus is not reproduced here. We are led to the identification of the zeroth order in p, r 



e a+da (P,R) = e Q (P,R)+ [d a U a (P,R)U+ (P,K) ,e a (P,K)]da 

i {A* l V Rl e a (P,R) + V Rl e a (P,R)^' +A%V Pl e a (P,R) + V Pi e Q (P,R)^} da 

{[e a (P,TL),A*>]A%da- [e a (P,R) t A%] A* - [AsymA™ ,s a {P,TL)]}da 



i 

2 

ida 



y {A S ym{V Pl V Rl e a (P,R)}-U a A S ym{V Pl V Rl Ho(P,R)}U+} (12) 

where we have defined the Berry phases [f| (note the absence of a in factor compared to the usual definition) 

A%> = iU a (P,R)V Pl U+(P,R) 
A% = -iU a (P,R)V Rl U+(P,R) 

and the "second order" Berry phase 

K" Rl = [Vfl fc V Pi [/ Q (P,R)] U+ (P,R) 

We have moreover introduced the linear operation Asym [5|. It acts on a symmetrical function in P and R in the 
following way : 



n 



Asym | -A (R) B (P) + -B (P) A (R)| = [fl (P) , A (R)] (13) 

the functions A (R) and B (P) beeing typically monomials in R and P arising in the series expansions of the physical 
quantities. 
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Note that the Asym operator especially in the term U+Asym {Vp, V Rl H (P,R)} [/+, is not particularly elegant. 
It could in fact be rewritten as an infinite series of differential operators, but the gain does not seem to be clear at 
this point. Moreover, in our examples, our choice will appear to be the most convenient. 

Ultimately, we arrive at the following differential equation 

^U Q (P,R) = [9 Q C/ Q (P,R)t/+(P,R),e a (P,R)] 

+ {^a ! V fli ea (P,R) + V Rl e a (P,R)^ +A%V Pl E a (P,R) + V Pl e a (P,R)^| 

+ l -{[e a (P,R),A*]A%- [e a (P,R) ,A%] A R ' - [AsymA p a - R \e a (P,R)]} 

+ l - {Asym{V P y Rl e a (P,R)} - U a Asym {V Pi V Rl H (P,R)} U+} (14) 

To complete the diagonalization process, we have to couple this equation with the evolution as a function of a of the 
transformation matrix U a (P,R). To find this last one, let us combine Eqs. [6l and [Til to get 

= U a (P, R)(7+(P, R)-U a (P + p, R + r)U+(P + p, R + v)+da [d a U a (P, R)C7+(P, R)+U a (P, R)C/+(P, R)] 

This latter expression can be, similarly to the energy, expanded as 

= Q (7 Q (P,R)[/+(P,R)+[/ Q (P,R)a Q C/+(P,R) - % -Asym (a™ + (A P ' R ') + ) - % - [A R ',A P '] (15) 

With these two equations Eqs. [JJ] and [J5] at hand, our diagonalization process can be performed. Actually, since 
all quantities are matrix valued and since e a (P, R) is by definition a diagonal matrix, we can separate the energy 
equation Eq. [14] in a diagonal and a non diagonal part such that we are led to the following two equations 

-^-£ a (P,R) = P+[R.H.S. of Eq. HI (16) 
da 

= 7>_[R.H.S. of Eq. [14] (17) 

where V+ [...] and V- [...] denote the projection on the diagonal and non-diagonal part respectively. These two equations 
are supplemented by the differential unitarity condition 

= a a C/ Q (P,R)C/+(P,R)+C/ Q (P,R)a Q C/+(P,R) - U S ym(A p ' R ' + «"*') + ) - \ [A Rl , A p ' ] (18) 

We claim that those three equations Eqs ll6l[T8l allow to determine recursively in powers of a the energy of the 
quantum system in question. Actually, the integration over a of Eq. 1161 gives e a (P, R) at order n in a when knowing 
all quantities at order n — 1. By the same token, Eqs. [T7] and [T51 (whose meaning is that U a (P,R) is unitary at 
each order in a) involve d a U a (P, R), and allow to recover U a (P, R) at order n by integration over a. Note however 
that these two equations allow only fixing partially U a (P, R), leaving m real parameters free with m the size of the 
matrix. This reflects a kind of gauge invariance (we will see a practical example of this in section four) . Actually, the 
diagonalization matrix is not unique, so that we are left with some choice of the parameters. For example one can 
impose the diagonal elements to be real (see below in the example). As a consequence, the diagonalization process is 
perfectly controlled order by order in the series expansion in a. 



III. THE SEMICLASSICAL APPROXIMATION 



In this section we consider the Hamiltonian diagonalization at the semiclassical level and the resulting equations of 
motion. Actually, the semiclassical approximation has recently found new important applications in particle and solid 
state physics. Notably, the equations of motion reveal a new contribution coming from the Berry curvature. This 
contribution, called the anomalous velocity, modifies profoundly the dynamics of the particles. For instance, the s pin 
Hall effect of electrons and holes in semiconductors Q , as well as the new discovered optical Hall effect [|[ 0] [HI [HI 
can be interpreted in this context. Similarly, the recent experimental discovery of the monopole in momentum can 
also be elegantly interpreted as the influence of the Berry curvature on the semiclassical dynamics of Bloch electrons 
00- 
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A. The semiclassical energy 



The consideration of Eq. [14] alone is sufficient to deduce the semiclassical diagonal Hamiltonian. Indeed, writing 
£a = £o + with e the diagonalized energy at the zero order, Eq. [T3]is solved by (putting a — K) 



e(P,R) = e (P,R) + fi. 



-A?V Rl e (P,R) 
ifi , 



X7 Rl e Q (P,R) A R '+A^'V Pl e (P,R) 



Vp ; e (P, R) A P ' } + l -V+ { [ Eo (P, R) , A&\ < ! - [e (P, R) , < ! ] < ! } 



where we have introduced the notations A R = V+ [A R ] and A^ = V+ [Aq] ■ 
This latter expression can also be written 



e(P,R) ~e(p,r) 



ih. 



An 



s(p,r),A Pl 



o(h 2 



e(p,r] 

where we have defined the projected dynamical operators (covariant coordinates and momentum operators) 

r = R + HAl 
p = V + hAl 



1? 



(19) 



(20) 



(21) 



with Af = i [U V-pU+] ,Ao=~i [UoV-rK] ,and A*' R = [V R V P [/ ] U+. 

The matrix £7o(P,R) is the diagonalization matrix for Hq when the operators are supposed to be commuting 
quantities, the diagonalized energy being eo (RR) • When P and R do not commute, the matrix Uq (P,R) does 
not diagonalize Hq anymore. In order to get the corrections to the energy at the semiclassical order due to the 
noncommutativity of P and R we have to compute e\ (P, R). Note that in Eq. [THl the symmetrization defined above 
(that is all the powers of P have been put half on the left and half of the right in each expression) is assumed. From 
the diagonal Hamiltonian, we can now derive the equations of motion for the covariant operators. 



B. The equations of motion 



Given the Hamiltonian derived in the previous subsection, the equations of motion can now be easily derived. As 
usual, the dynamics equations have to be considered, not for the usual position and momentum, but rather for the 
projected variables r and p . Actually, these latter naturally appear in our diagonalization process at the H order. 
Let us remark, as now well known, that their components do not commute any more. Actually (removing the index 
a for clarity) 



[ n ,rA = ih 2 QY 3 = rt 2 (Vp.4, - Vp^rJ + h 2 [A Rj ,A Rt ] 

[ Pt , P3 ] = m 2 e^ = -in 2 (v Ri A Pi -v Rj Ap i )+h 2 [Ap i ,Ap j ] 



-ihSi. 



ih 2 (V Ri A Rj + V Pj A Pi ) + h 2 [A Pi ,A Rj ] 



(22) 



the @ij being the so called Berry curvatures. 

Using now our Hamiltonian yields directly to general equations of motion for r, p 



f = I[r, e (p,r)] + I 

i i 
P = j: [P,e(p,r)] + - 



ih 

r, y^+[[ £ (P,rM««Mfl - [e (p, r) , A Pl ] A Rl ] 



P,y^+ [[e (P,r),A Rl ]A Pl - [e (p, r) , A Pl ] A Rl ] 



(23) 



The commutators can be computed through the previous commutation rules between r and p. The last term in 
each equation represents a contribution of " magnetization" type and has the advantage to present this general form 
whatever the system initially considered. In the context of Bloch electrons in a magnetic field, it gives exactly the 
magnetization term revealed in [l[ (see @). For spinning particles in static gravitational fields, this term gives a 
coupling between the spin and the intrinsic angular momentum with magneto-torsion fields 14j 15(. It is interesting 
to note that both the Hamiltonian Eq. [20l and the equations of motion have already been derived in article [5] .Our 
method appears then to be an extension of [5|. However, although relying on this previous method, it seems to us 
more elegant, and offers to possibility to go beyond the semiclassical level. This will appear very important in the 
next section to compute the velocity of a spinning massless particle (photon) in an inhomogeneous medium. The 
second application is the determination of the equations of motion of a Bloch electron in an external constant electric 
field. In both situations we are going to consider the second order in h quantum effects. 
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IV. PHYSICAL APPLICATIONS 



A. Massless spinning particles in an isotropic inhomogeneous medium. 

We consider here the dynamics of a spinning massless particle propagating in a isotropic inhomogeneous medium 
of index n(R). We thus start with the Hamiltonian of a massless Dirac particle in an curved space of metric 
g y (R) — rT 1 (R) (5 U as in The case of higher spin (like the photon) can be treated in the same manner 

through the Bargman-Wigner equations as in [l4| . For this reason we limit ourself to the Dirac particles (with the 
convention c = 1) and write the Hamiltonian as 

H = - («.PF(R) + F(R)a.P) (24) 

with -F(R) = n^ 1 (R). Let us start to recover the first order diagonalization given in To exploit the previous 
general formula we write all quantities as a series expansion in a, thus e a = £q + ae\ + U a = Uq + aU\ + ■ ■■,A^ — 
Aq + aAi + ... and the same for all expressions. In these expressions, all terms are functions of P and R and are 
thus assumed to be symmetrized in the sense defined above. 



1. First order (semiclassical) diagonalization. 



We already know (see [14|]) that at the zeroth order in a, that is when R and P commute, the Hamiltonian 
diagonalization can be performed through the following Foldy-Wouthuysen unitary matrix 



U (P) 



(3a.V 



/2p2 



(25) 



(where j3 and a are the usual notations for the "Dirac matrices" ) such that in the diagonal representation the 
Hamiltonian becomes 



U H (P, R) U+ = e Q (P, R) = \ (/3F(R) VT^ + /?VP^(R)) 



(26) 



From the unitary matrix, we can deduce that at this order, only „4 R = i [C/Vp[/ + ] is non zero, thus .4^ = 
— ih [C/Vr?7 + ] = + O (a) ,and Aq 1 ' Ri = + O (a) . Obviously, when P and R do not commute, the matrix e Q (P, R) 
does not diagonalize Hq anymore and one has to go to higher order in the expansion in a. 

Let thus consider the diagonalization at order a. The differential equation Eq. [I4]at the zeroth order (at this order, 
we consider only in the right hand side of this equation the terms of order a ) reduces to 



jU Q (P,R) = i [< ! V fll e (P,R) + V fi;£o (P,R)<'] + [t/i(P,R)[/ + (P),e (P,R) 



and can be divided in two parts 



and 



= -V. 



(P,R) - ip + [A*'V Rl e (P,R) + V Bi e (P,R)<' 



A*' V Rl e (P, R) + V^eo (P, R) A$'} + [U x (P, R) U+ (P) , e (P, R)] 



(27) 



(28) 



(29) 



(the last term being always anti-diagonal at this order). By integrating Eq. [55] with respect to a from o to H the first 
equation, we get the energy operator at the first order ina^li which reads 

e a (P, R) = £ (P, R) + ^V^eo (P, R) + ^V fl ,e (P, R) < ! 

The energy e a can be recombined under the usual form [j| and written in term of the covariant coordinate operator 
r+ as 

e Q (P,r) = !(F(r)\/p^+\/p^F(r)) (30) 
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with 

r = R+hA* = R+ifi P * 2 S (31) 

where means that we need only the diagonal part of „4^. The matrix S = l(g)crisa4x4 matrix. 

From the last two equations we can deduce the equations of motion in the semiclassical approximation. Indeed, for 

a particle of positive energy only e (P, r) = | (f (r) a/P^ + \fP 2 F (r)j with now r = R+i^p^ a 2 x 2 matrix (the 

spin matrix is S = ha/2), the usual relations r = —i [r, e(P,r)] and P = — i [P,e(P,r)], lead to the by now well 
known equations of motion 

f = Vp£ + fiPx8" 

P = V r e (32) 

where [r^r?] = ih 2 Q r ^ = ih 2 eijk&l r — ~~ ihzijk^^pz, and A = S.P/P the helicity. These latter equations are very 
important because they induce the optical Hall effect or spin Hall effect of light (for the photon one just has to replace 
the Pauli matrices a by the spin one matrices) as a consequence of the anomalous velocity term hPxQ rr . This 
effect has been discussed in several circumstances [l(| [H| and was considered in a quantum mechanical context in 
000. From Eq[32] we see that the velocity components are given by the relation 

i _ 1 ( c P i P^^_\ _A c_ k dhm 

V ~ 2 \n(r) P + P n(r) ) + P 2 n(rf ljk da* [ ' 

Note that the symmetrization for the last term was unnecessary due to the presence of the h. Symmetrizing would 
just add H 2 corrections. 

From the last relation, we deduce the velocity 

l + ^((Vlnn) 2 -4(P-Vlnn) 2 ^ ' (34) 



n(r) V P 2 V P 2 

The velocity can be rewritten v = (1 + -p^ (Vlnn) sin 2 i9J with i9 the angle between P and Vlnn. We 

therefore find the important result that the light velocity in an isotropic inhomogeneous medium does not have the 
well known expression v(r) — c/n(r) but has rather quantum corrections due to the spin-orbit couplings when Pand 
Vlnn are not parallel. However, at this level of approximation, we can not yet conclude regarding the velocity since 
the quantum corrections are of order h 2 which is beyond the first order approximation. To get a clear answer one has 
to compute the Hamiltonian diagonalization at the second order in h. 



2. Second order diagonalization 



To perform the diagonalization at order a 2 we need to know the transformation matrix U a at the first order in a, 
i.e., the matrix C/i(P,R) (i.e. the Foldy Fouthuysen transformation when R and P do not commute). The equation 
in Eq. [21] allows to find easily the transformation matrix U a at this order. Indeed, given that at zeroth order 



V- 







A* 1 V Rl e (P, R) + Vh, e (P, R) 1 

as a consequence of 5/3 + /35 = 0, we can write 

[[/ 1 (P,R)[/ +(P),eo(P,R)] =0 
Likewise, the differential unitarity condition Eq. [TS]at the same zero order becomes 

iMP,R)[/ + (P)+;7o(P)t/|(P,R) =° 



(35) 



(36) 



(37) 



We claim that U a (P, R), and consequently the Berry phases, are unchanged with respect to the zeroth order. That 
is the matrix Eq. [55] also diagonalizes Ho at the first order in a. Indeed, we can easily see with the help of Eqs. I3AB71 
that the following transformation (l — aUi(P, R)Uq(P)) C/ q (P,R) which at order a is equal to Uq(P) diagonalizes 



9 



Hq, e.g. UqHqUo = £o (P,R) + asi (P,R) (this is an example of the gauge transformation in the diagonalization 
matrix we mentioned previously). 

As a consequence, the Berry phases do not get any contribution at the first order. Therefore we have = 
-iV+ [UW R U+] =0 + O (a 2 ) implying A%' Rl = + O (a 2 ) and also 



AK = i hV + [UV P U+]= i ^ + 0(a 2 ) 



Now, to find the second order Hamiltonian we write the differential equation Eq. rjj]at the first order. 

-e a (P,R) = \v+ [A Rl V Rl e a (P,R) + V Pi £ Q (P,R)A Rl ] + l -Asym {V p V P ,£ Q (P,R)} 



(38) 



(39) 



since trivially Asym {V p W R[ H (P, R)} = 0. To show how to solve this equation, recall that at the first order 
e a (P, R) is given by 



e a (P, R) = £o (P, R) + f A* V Rl e a (P, R) + |v Pi £ a (P, R) A* 



(40) 



(recall that all our expressions are supposed symmetrized in P and R as explained in the beginning of the previous 
section), so that Eq. I3"9l becomes 

-^e a (P, R) = \A%V Rl ev (P, R) + iv fll e (P, R) A R > + ~A R « A* V Rk V Rl e (P, R) 

+f^ h V Rfc V Ri e (P,R)A*< + |v Rfe V Ri£o (P,R) A* A* + l -Asym{V P V Rl e (P,R)} 
and as a result after integration, 

2 

e a (P,R) - £o (P,R) - - A*V Rl e (P,R) + |v Pl £ (P,R) A R > + —A R "A R ^ Rk V Rl e (P,R) 

2 2 

" A R *V Rk V Rl e (P,R)A R ' + ^-V Rk V Rl e (P,R)A R ><A R > 



4 

'2J0 



daAsym{V Pi V fli e (P,R)} 



Given our symmetrization convention, this last expression e a (P, R) can be shown to be equal to £o (P, r). 

Given this last result, we can solve our differential equation at the second order approximation. Actually, recalling 

that £ (P, R) = | {j3F{R) VP 2 + /?VP~2F(R)) , one can first state that : 



A S ym{V Pi V Pi £ (P,R)} = VP2,F(r) 



so that 



i 

2 J 



da Asym {V p V Pl £ (P, R)} =--H [VP 2 , F (r) 



Now, by replacing a by h and projecting on the positive energy subspace, our solution is :. 



£(P,r) = i(V(r) 



1 



/pa + (r) ) - -h \yV 2 ,F (r) 



which can be also written 



£ (P, r) = i (f (r) VP 2 + Vp^f ( 



1 H 2 

)'! ! - t— PVf(r) 



4 VP2 



(the symmetrization of the last term is not necessary due to the presence of the ti 2 .) 

Now, concerning the equation of motion at the second order, their form is only formally identical to Eq. [32] 
Actually, the energy has been changed by getting quantum corrections of order h 2 . Moreover, the term HP x Q rr in 
EqJ32]has to be understood as beeing symetrized between P and r. 
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Let us conclude this section by noting from the equations of motion that the speed of light get also quantum 
correction of order f? . Indeed, one has now a correction to Eq. [33] of order H 2 (remember that A is of order H) such 
that the velocity components become 

I / c pi P l c \ A f ,,d\nn c Sinn c A ti 2 c / 1 „ , 



2 Vn(r) P ' P n(r) J ' 2P^ m \ P dx* n(r) + dx< n(r) P J + 4 n(r) ^p 9 ' 1 ™ P 3 ^ lnn 

Interestingly the corrective term does not give any contribution to the velocity at this order of approximation so that 
the velocity is indeed given by Eq. [3J] at the order ti 2 

\2 / I \ \ 1/2 

A / , ,2 1 , . \2 



1 + -= (V Inn) (P. V In „)* + 0(fi d ) 



n(r) V. -P 2 V p2 

It is worth mentioning that in the context of photons propagating in a smoothly inhomogeneous medium we are 
apparently led to a contradiction with known results. Indeed, it is known from optics that the second order corrections 
in a smoothly inhomogeneous medium are responsible for the linear birefringence, i.e., the equations of motion depend 
on the linear polarization of the particle [18j and not only on helicity as found here. One possible explanation of this 
discrepancy might be that here the photon in an inhomogeneous medium is treated as a quantum particle moving 
in a static gravitational field, so that the medium index is equivalent to a metric. This approach is perhaps too 
restrictive since by doing so, we take only into account the geometry and neglect all electrodynamical processes which 
are relevant for the interpretation of the linear birefringence [l8j |. In other words, an inhomogeneous medium might 
be considered to be equivalent to a gravitational field at the semiclassical order only, since the electrodynamics plays 
a role at the second order. If true, our approach beyond the first order is strictly speaking only valid for spinning 
particles in a static gravitational field. 



B. The electron in a time dependent electric field 



As a second example, we consider a spinless non-relativistic electron in a periodic potential submitted to an external 
time dependent electric field. We limit ourself to this simple case because a second order computation can easily be 
performed as shown below, and leave more complicated cases with a position dependent electric field and a magnetic 
field to future works. Because of the non-trivial vacuum configuration, that is the band structure of the energy 
spectrum in solid state physics this example fits in our formalism similarly to the Dirac equation. Important new 
results where recently found in the context of Bloch electron. Indeed, the semi-classical equations of motion where 
found to be modified by Berry phase terms [l[ p| 0] [l9| , changing profoundly the properties of electron transport in a 
solid 20]. It is thus interesting to go beyond the semi-classical approximation to check wether new effects might be 
revealed. 



1. First order 



Since the initial diagonalization matrix when a = is not know, the results obtained are more formal. Nevertheless, 
let us start with the 0th order formal diagonalization. Starting with the Hamiltonian of an electron in a periodic 
potential plus an external perturbation 

ff (P, R) = ^+^(R)+«(R) (41) 

with v(R) = E.R and E(£) a time dependent electric field, the zeroth order diagonalization can be done formally 
with a K dependent unitary matrix U a (K), where K is the pseudo- momentum of the Bloch bands @ [|[. Therefore, 
only the position acquires a Berry phase, i.e., (K) = i [UWpU + ] + and A& = —i [UWpU + ] + =0 + O (a) implying 
A%' R ' = + O (a) . 

As a consequence, the diagonalized Hamiltonian e Q (K, R) at order a, projected on a particular band of index n is 
given by : 

£ ai „(K,r n ) =e n (K)+«(r n ) (42) 

where e n (K) is the energy of the unperturbed nth band (i.e. when w(R) =0) and r„ the covariant position operator 
defined as the projection of R + aA R on the nth band: 



r„ = R + o4£ Q (K) 



(43) 
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where A^ n is the projection of the matrix A^ on the nth band. The equations of motion are easily deduced to be 

t„ = V K £ Q ,n + ^Kxe; r 

K = V r e Q ,„ = E (44) 

with the Berry curvature defined as [r n ,j, r n ,j] = ©JTij = ^iAn, a ~ 9jA^ a . The second term in Eq. 25] is called 
the anomalous velocity and it is was shown todrive important new effects, such as the discovery of a monopole in 
momentum contribution to the conductivity of Bloch electrons in particular materials (l6| . It is then interesting to 
check if this equation is changed beyond the semiclassical level. 

2. Second Order and more 

The computation of the energy to the second order needs the value of the transformation matrix at the first order 
in a. Consider therefore the Oth order differential equations for the transformation matrix Eq. [17] 

= [d a U a (K, R) U+ (K, R) , e a (K, R)] ~ + ~ [A*>V Rl e a (K, R) + V Rl s a (K, R) (45) 

and Eq. Ql 

= d a U a (K, R)£/+(K, B.)+U a (K, R)d a U+(K, R) (46) 

We see that nothing ensures that : [A^ 1 V_r ; £ q (K, R) + V^Eq, (K, R) A^ 1 ] = 0. To solve for the transformation 
matrix, we proceed as for the photon and decompose at the first order E/ a (K,R) =t/o(K, R)+aJ7i(K, R), so that, 
taking into account the specific form of t?(R) our equations Eqs. l45l and l46l reduce to 

[E7iCK,R)E/ + (K J R),e a (K,R)] = -[A^,Ei]~ 
[/i(K,R)C/ + (K,R)+[/i(K,R)f/+(K,R) = 

Introducing the band indices, one easily sees that these two equations reduce to the 2n 2 — n following independent 
equations for 2n 2 real variables (n, the number of bands can of course be infinite; this does not invalidate the argument) 

[l7i(K,R)Ctf(K,R)] mn (e a (K,R)„ - e a (K,R)J = - [A*,®]^ , for m<n 
[CT 1 (K J R)t^(K ) R)] mn +[l7i(K J R)C^(K J R)] Tnn = 

Once again, we are left with a gauge choice on the diagonal. We can simply impose as a condition that the diagonal 
elements are null. Actually, our equation imposes that [f/i(K, H,)Uq(K, R)] is an imaginary number. As a 
consequence, similarly to the photon, composing C/ Q (K, R) on the left with the unitary matrix X defined by 

X mn = 1 - aS mn [tfi(K, R)U+(K, R)] m (47) 

yields a diagonalization matrix with zero on the diagonal. 

Having chosen the gauge, the equations for £/i(K, R)C/q (K, R) can be solved formally, through the Cramer rules. 
Fortunately, given our very simple case of potential, we get 

e a (K, R)„ - e a (K, R) m = s a (K)„ - e a (K) m (48) 

which does not depend on R. As a consequence, t/i(K,R)[/ + (K,R) and thus [7 X (K,R) is only a function of K. 
We can thus write 

tf a (K,R)=Z7 (K)+al7i(K) (49) 

As a consequence at the first order we will have again.: 

A^t[U a W K U+] + + 0(a 2 ) (50) 

and for the other Berry curvatures we have A^ =0 + O (a 2 ) , and A^ l ' Rl =0 + O (a 2 ) . Let us just mention that 
the Berry phase in R is now different from its the zeroth order value. It has now a term proportional to a which is 
included in its definition. 
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Consider now the differential equation for the energy at the first order in a. A simplification arise, since here, at 
the first order 

Asym{V P ^ Rl e a (K,R)} - U a Asym { V Pl V Rl H (P,R)}C/+ = (51) 

actually K and R are not mixed in e a or Hq. As a consequence, and given our previous results, the projection of our 
equation on the diagonal part leads to 

^-e a (K,R) = 1 [A B a l V Rl e a (K, R) + V Rl e a (K,R)^] + (52) 

Exactly as in the photon case, this equation has the second order solution after projection on the nth band 

e a . n (K,rJ = e„(K) + v (r„) (53) 

where now r„ = R+aA^ n - Note that aA B n is actually of order a 2 . Therefore the energy has formally the same 
expression than in the first order case. The only difference stems from the potential energy because the position 
operator get a a 2 contribution through the presence of the Berry connection. 

As a consequence, the equations of motion are formally the same as for the first order (removing the band index), 
i.e., 

1 • 

r„ = V K £a,« + ^ K xQ r n r 

m 

K = Vr n £ Q ,n = E 

The only difference with respect to the semiclassical equations of motion Eq[44] is that the Berry connection gets a 
contribution of order h 2 but formally the equations of motion are unchanged. 

To conclude this section let us state the following important result. One can easily check that all our construction 
can be generalized to all orders in h. As a consequence, the Hamiltonian form for e a , n (K, r n ) is valid to all orders, and 
the equation of motions are in fact quantum operator equations. The only difference appearing when diagonalizing 
recursively, consists in new contributions to the Berry connection. 

This is a important new result, because these equations were considered until now has semiclassical equation of 
motion with the usual restrictions imposed in this case, ft would be interesting to check if this result is generalizable to 
spatially variable electric field and in the presence of a magnetic field, but then the equations get quickly cumbersome. 

V. CONCLUSION 

Some recent applications of semi classical methods to several branches of Physics, such as spintronics or solid state 
physics have shown the relevance of Berry Phases contributions to the dynamics of a system, leading notably to the 
discovery of the intrinsic spin Hall effect. However, these progresses called for a rigorous Hamiltonian treatment that 
would allow for deriving naturally the role of the Berry phase in theses systems and to go beyond the semiclassical 
level. In this paper, we presented a new diagonalization method for a generic matrix valued Hamiltonian, based 
on a formal expansion in power of h. A differential equation connecting two diagonalization processes for two very 
close values of h, considered as a running parameter, was derived. The integration of this differential equation allows 
the recursive determination of the series expansion of the diagonalized Hamiltonian in powers of h. This approach, 
which results in effective Hamiltonians with Berry phase corrections, goes beyond previous works on the semiclassical 
diagonalization of quantum Hamiltonians. The resulting generic equations of motion are also corrected by Berry 
phase terms of higher order in h. As physical applications, we considered a spinning massless particle propagating 
in a isotropic inhomogeneous medium and showed that both the energy and the velocity get quantum corrections 
of order h 2 . We also derived formally to all order in h the energy spectrum and the equations of motion of Bloch 
electrons in an external constant electric field, showing that the equations of motion are actually quantum operator 
equations (valid to each order in H). It would be interesting to check if this result persists for variable electric fields 
and in the presence of magnetic fields. Our approach is a general one and will be applied in the future to other 
condensed matter systems as well as in particle physics. 

Acknowledgment. The authors acknowledge a fruitful correspondence with K. Bliokh about the linear birefrin- 
gence phenomenon discussed in the section IV. 2., as well as for having brought to our knowledge reference [18| 
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